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1. Introduction . Active galactic nuclei (AGNs) and quasars (QSOs) appear 
to emit roughly equal energy per decade from radio to y ~ ray energies 
(e.g. Ramaty and Lingenfelter 1982). This argues strongly for a non- 
thermal radiation mechanism (see Rees 1984). In addition, statistical 
studies have indicated that the spectra of these objects in the IR-UV and 
2-50 keV X-ray band, can be fitted very well with power laws of specific 
indices. These spectral indices do not seem to depend on the luminosity 
or morphology of the objects (Rothschild et al. 1983; Malkan 1984), and 
any theory should account for them in a basic and model independent way. 

If shocks accelerate relativistic protons via the first-order Fermi 
mechanism (e.g. Axford 1981), the radiating electrons can be produced as 
secondaries throughout the source by proton-proton (p-p) collisions and 
pion decay, thus eliminating Compton losses (Protheroe and Kazanas 1983). 
As shown by Kazanas (1984), if relativistic electrons are injected at 
high energies, e -e pair production results in a steady state electron 
distribution that is very similar to that observed in AGNs, independent 
of the details of injection and the dynamics of the source. The 
conditions required by this mechanism are met in the shock model of 
Eichler (1984) and Ellison and Eichler (1984) which allows the self- 
consistent calculation of the shock acceleration efficiency. 

2. Model . We assume that a black hole of mass, Mg = M/(10 9 M), accretes 

gas of density, n, and temperature, T. If the dissipative tfme scale of 
the accreting material, t , is longer than the free fall time, x ff , a 
steady state radiative shock will develop at x ■> R/R , where R 
= 2GM/c 2 is the Schwarzschild radius. Inflowing thermal particles with 
free fall velocity, u, will be shock accelerated and relativistic protons 
will be produced. These protons will undergo inelastic nuclear collisions 
to produce pions, i.e., p + p p + p + + 7 r°. The pion decay will 

produce yrays, neutrinos, and relativistic electrons which will 
subsequently produce the observed radiation by synchrotron and inverse 
Compton emission. Since ~ 50% of the energy flux that goes into 
relativistic particles is lost to neutrinos (Eichler 1979) and since the 
shock must form beyond the Schwarzschild radius, the maximum efficiency 
for this model is considerably less than 1. The maximum luminosity is 
determined by the ability of the shock accelerated particles to provide 
sufficient pressure to support the shock against the accreting material. 
This is determined by the shock dynamics and occurs at ~10% of the 
Eddington luminosity. 

A fraction, Q, of the energy flux crossing the shock is converted into 
relativistic particles and eventually into radiation. Therefore, 

L tot " 4 * R i * i n r P u i • q. a) 

where L to(; is the total radiated energy including neutrino losses and the 
subscript (1) indicates values at the shock. Due to neutrino losses, the 
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observed luminosity, L, will be L ~ 0.5 L tot . Alternatively, if <E re] > is 
the average energy density in relativistic particles, the total 
luminosity can be expressed as, 

L tot ■ V- (2) 

3/2 2 —1 

where x - l/(<n,^oc) - 2,2 10 8 xj Mg m sec is the p— p interaction 
time, cP%.s the p-p cross section, xj is the shock position in Schwarz- 
schild radii, <n ± > is the average density inside the shock, and m is 
the accretion rate in M Q /yr. Equations (1) and (2) determine t>he shock 
position. Assuming mass conservation, i.e., n = 1.2 10 9 x Mg m/cm 3 , 
and using <E rel > = 3P rel , where P rel is the relativistic particle 
pressure, we have 



(3) 


where n = P rel /(nim uj 2 ). Eq. (3) indicates that m/M 9 is determined 
solely by the shock Efficiency and cannot have an arbitrary value. The 
efficiencies Q and n (which to good approximation depend only on the 
shock Mach number) can be calculated self-consistently using the non- 
linear, steady state solution of Eichler (1984) and Ellison and Eichler 
(1984). This assumes that a magnetic field is present for the formation 
of the collisionless shock and the elastic scattering of the accelerated 
particles. We assume, however, that the magnetic field is not dynamically 
important and that the AlfvSn (or phase) velocity is small compared to 
the shock velocity. We also require that the mean free path of the 
highest energy particles be no smaller than R g so they avoid being 
swallowed by the black hole and hence can provide the pressure to support 
the standing shock. 

Since both Q and p approach 1 for high Mach numbers, there is a 
maximum value above which no steady state shock solutions of the type 
described here exist. If m/M 9 is greater than this limit, the p-p 
collisions behind the shock occur so rapidly that there is not sufficient 
relativistic particle pressure to support the radiative shock. 

Using the above relations, the observed non-thermal luminosity , L, 
becomes 


L « 1.5 x 10 97 Q 2 M 9 /(n x i) erg/s. (4) 

, Defining, F = L/L E , where L £ “ 1.3 x lO^ 7 M 9 erg/s, eq. (4) reads 

F » L/L e « 1.2 Q 2 /(n xj). (5) 

* This indicates that the sources should emit at a given fraction of the 
Eddington luminosity (L E ), which depends only on x x , Q, and p. 

The value of x x can be expressed in terms of the Mach number, Mj, and 
the upstream temperature as, 

M^ = 3m p u^/[5k(T e + T p )] « 6.5 x 10 4 /(x L T g ), (6) 

where T 8 is the upstream temperature in units of 10 8 K and is expected to 
be determined by the balance between X-ray Compton heating and cooling 
and has been shown to be essentially constant and ~ 10 8 K (e.g., Krolik 
et al. 1981). Consequently, F, depends on only one parameter, Mi or Xj. 



1 30 


0G 8.1-7 


3. Results and Conclusions . By combining eqs. (5) and (6), we can obtain 
the Eddington efficiency as a function of the shock size or the accretion 
rate. These relationships are shown in Figs, la, and b. A range of 
upstream temperatures is shown. 

The Eddington efficiency 
drops precipitously with in- I 

creasing x^ due to the fact that 
Q approaches zero as M 1 drops |Q -i 

below ~ 4. One therefore expects 
that there would be an effective 
upper limit for xj (lower limit 10 2 

for Mi) above which the sources 
are too weak to be seen. At the io -3 

other extreme, there is a max- * 
imum value for F since Xi can 
not be arbitrarily small. We 10 

have chosen x^ = 5 as a 

reasonable minimum value of the !0' s 

shock radius. This corresponds 
to an efficiency, F =0.2. 10 - 6 

We can also express the lumi- 
nosity as a function of the 
black hole mass with xj as a x, m/M 9 

parameter. Our results can now Fig. 1. L/Lg vs. shock position 

be directly compared to the x*, and accretion rate, m p 

correlation between L and M Lg/c^ — 2.3 Mg solar masses/yr. 

(Fig. 2) found recently by 

Wandel and Yahil 1985) (W-Y). These authors, by attributing the widths of 
the Hg line to dynamical motions, and using photoionization arguments to 
determine the distance of the line emitting clouds, showed that the 
luminosity of ~ 70 QSOs and Seyfert 1 galaxies is proportional to their 
mass. In comparing our model, which yields the entire non-therraal 
luminosity, a bolometric correction raising the points of W-Y by 5-10 
should be used. With such a correction, the shock model with 5 £ xj, £ 150 
reproduces the observations extremely well. The observations^ indicate 
that even though the non-thermal radiation from QSOs and AGNs is 
generally well below the Eddington limit, it is still highly correlated 
and proportional to the mass. This proportionality is an integral feature 
of our shock model and results from the fact that m « M g (eq. 3). The 
proportionality constant depends mainly on the p-p strong interaction 
cross section and results in luminosities well below the Eddington limit 
in accordance with observations. Also, the fact that no objects in the W- 
Y compilation are observed to emit much below 0.01 L E might have an 
explanation in the sharp decrease in F with m/M g as indicated in Fig. lb. 
We consider the straightforward interpretation of this correlation by our 
model as an indication of the correctness of its basic premises. 

The galaxy NGC 4151 affords an additional test of our model. The 
considerations of W-Y provide an estimate of the mass, M = 3 x 10 7 M . 
This mass along with the observed absolute luminosity, L = IO 43 erg/s, 
determine, F = 2.5 x 10 -J . This is shown by the horizontal dashed line of 
Fig. la labeled "dynamics”. X-ray variability has been established on 
time scales of At = 12 hours (Mushotzky et al. 1978), indicating a size, 
R = cAt = 1.3 x 10 15 i cm. Using the determined mass, this translates into 
a shock radius, x^ = 140, and is shown in Fig. la as "time variability". 
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Furthermore, measurements of the X-ray spectrum of NGC 4151 (Mushotzky et 
al* 1978) provide a constraint on the compactness, L/R, within the 
framework of the non-thermal, e + -e~ feedback model of AGN spectra 
(Kazanas and Protheroe 1983; Kazanas 1984); i.e., L/R £ 10 28 erg/(s-cm). 
This condition becomes, F £ 2.3 x 10 x^, and is indicated in Fig. la by 

the line labeled "spectral information". We find the agreement between 
these independent constraints and our model to be quite remarkable. 

In conclusion, we have 
outlined a model for AGNs and 
QSOs which is complete in the 47 1 
sense that it simultaneously 
addresses the source dynamics 
(i.e. the conversion of accretion 
energy into radiation) and 
spectrum (from radio to Y“*ays)< 

As shown in Kazanas (1984), the Jfc 
e + ~e w feedback decouples the ^5 
spectral signatures from the >f 
dynamics, provided that most of to 
the power is injected at energies 2 44 
high enough that the resulting § 
photons are absorbed within the ^ 
source to produce e + -e“ pairs. § 43 
The present model does precisely 
that, through shock acceleration 
(Ellison and Eichler 1984). 

Furthermore, it provides a 
straightforward physical argument 
as to why such shocks should 
occur (Xpp > Tff) and also gives 
estimates of the Eddington 
efficiencies, F, and sizes, x\ 9 
of these sources. Finally, in its 
most direct confrontation with 
observations (Fig. 2), our model provides a natural explanation for the 
L « M correlation. In addition, it also provides the normalization (it is 
directly related to x p ) and its value is in remarkable agreement with 
the data of W-Y. 
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2. The lines are from the 
shock model (using Tg » 1 ) for 

shocks which form at various x^* 
Points are from Wandel and Yahil. 
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